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Gautam C SethisQ and Abhijit Seij^ 

Institute for Plasma Research, Bhat, Gandhinagar 382 4-28, India 

(Dated: December 13, 2011) 

We investigate the existence and stability of traveling wave solutions in a continuum field of 
non-locally coupled identical phase oscillators with distance-dependent propagation delays. A com- 
prehensive stability diagram in the parametric space of the system is presented that shows a rich 
structure of multi-stable regions and illuminates the relative influences of time delay, the non-locality 
parameter and the intrinsic oscillator frequency on the dynamics of these states. A decrease in the 
intrinsic oscillator frequency leads to a break-up of the stability domains of the traveling waves into 
disconnected regions in the parametric space. These regions exhibit a tongue structure for high con- 
nectivity whereas they submerge into the stable region of the synchronous state for low connectivity. 
A novel finding is the existence of forbidden regions in the parametric space where no phase-locked 
solutions are possible. We also discover a new class of non-stationary breather states for this model 
system that are characterized by periodic oscillations of the complex order parameter. 

PACS numbers: 05.45.Ra, 05.45.Xt, 89.75.-k 



I. INTRODUCTION 



The classic Kuramoto model |Ij , consisting of a ring of 
phase coupled oscillators, has for long served as a very 
useful and convenient paradigm for investigating collec- 
tive phenomena in a variety of physical, chemical and 
biological systems [2 - 5j . The original model employing a 
simple global (all-to-all) and constant uniform coupling 
between the oscillators, has subsequently been extended 
and generalized in various ways to improve its applica- 
bility as well as to explore a wider range of dynamics [Bl- 
[15] . One such generalized version that provides a closer 
representation of many real life situations, employs non- 
local coupling to account for diminution in the coupling 
strength as a function of distance and distance depen- 
dent time-delayed coupling to account for finite propaga- 
tion velocities of information signals exchanged between 
oscillators [TS] . Such a model displays a much richer vari- 
ety of collective excitations including synchronous states 
[15j and exotic chimera states |14j . Both non- locality and 
time delay play important roles in influencing the equi- 
librium and stability properties of these collective states 
and their combined presence often introduces novel ef- 
fects. To the best of our knowledge, an analysis of the 
equilibria and stability of traveling wave states for the 
generalized model has not been done so far, excepting a 
few past studies that have employed simplified models or 
looked at special limits. Crook et al. T] considered a con- 
tinuum system of coupled identical oscillators with a spa- 
tially decaying interaction kernel and modeled the space 
dependent time delay contribution through an effective 
phase shift term in the interaction. They also took the 
size of the system to be infinite for mathematical conve- 



* Electronic address: g;autam@i pr.res.ln | 
tElectronic address: abhijitffipr.res.in 



nience but thereby effectively reduced the nature of the 
mutual coupling to a "local" one (since the interaction 
length is always much smaller than the system length). 
Zanette [5] studied another simplified version of the gen- 
eralized system where he adopted a distance independent 
(global) coupling between the oscillators positioned along 
a ring but introduced a distance dependent time delay in 
the interaction. He obtained numerical results on the 
stability of synchronous and propagating traveling waves 
and also some analytic results in the limit of small delay. 
Ko & Ermentrout [13] recently investigated the effects of 
distance dependent delays in sparsely connected oscilla- 
tor systems and found that a small fraction of connections 
with time delay can destabilize the synchronous states. 
Since traveling wave states are basic ingredients for pat- 
tern formation in many real life systems such as arrays 
of Josephson junctions [TB], chemical oscillators pUlIlT) . 
neural networks producing snail shell patterns and ocular 
dominance stripes |18j etc. a determination of their sta- 
bility properties within the framework of the generalized 
model is of fundamental importance. 

In this paper we present a systematic and comprehen- 
sive stability diagram of traveling wave states adopting 
a suite of numerical approaches that has proved success- 
ful in the past in delineating the stability of synchronous 
states of this fully generalized model [TS] . The stability 
diagram shows a rich structure of multi-stable regions in 
the parametric space of time delay and non-locality of 
the coupling with the traveling waves losing their stabil- 
ity across the marginal curves through a Hopf bifurca- 
tion. Interspersed among the stable regions there exist 
forbidden regions (whose size and shape show a strong 
dependence on the intrinsic frequency of the oscillator) 
in the phase diagram where no stable phase-locked solu- 
tions are possible. Our numerical investigations also re- 
veal the existence of a new class of non-stationary states 
for this system that are characterized by periodic oscilla- 
tions of the complex order parameter and for which the 



phase and its temporal derivative have a non-hnear de- 
pendence on space. Finally in support of our detailed 
numerical results we also propose an analytic heuristic 
necessary condition for stability that provides an upper 
as well as a lower stability limit for the traveling wave 
solutions. 



II. MODEL SYSTEM AND ITS TRAVELING 
WAVE STATES 

Our dimensionless model equation representing the dy- 
namics of a continuum of coupled identical phase oscilla- 
tors located on a ring, x £ [—L,L] is given by [l5] . 



where k is a dimensionless quantity denoting the inverse 
of the interaction scale length and is a measure of the 
non-locality of the coupling. 

Before any further deliberations it is convenient to de- 
fine a mean delay parameter by 



— wr — = w- / G{z) sin [4i{x, t) ~ (t){x - z,t - \z\t„i)] dz. 

(1) 
where (^(x, t) G [0, 27r) is the phase of the oscillator lo- 
cated at X and at time t and whose intrinsic oscillation 
frequency is w > 0. In ([I]) the space co-ordinates x and z 
have been normalized by L and time and frequencies have 
been made dimensionless by the prescription t = Kt, 
u! = Ld/K where K is the strength of the coupling be- 
tween the oscillators. The quantity Tm = 1/v denotes the 
maximum time delay in the system with v representing 
the signal propagation speed. The normalized function 
G : [— 1, 1] — >■ M describes the coupling kernel and is an 
even function of the form, 



G{z) 



2(1.0 -e-«) 



-K.\z\ 



(2) 



G{z)i 



I dz 



(3) 



which weights the individual delays with the correspond- 
ing connection weights and for the exponential connec- 
tivity given by Eq. ([2]), becomes 



= T„i/C 



(4) 



where c, 



_ ^(6"-!) 



We look for phase-locked solutions of Eq. ([I]) of the 
form 



4'n,kix, t) = Qt + nkx + (j)^. 



(5) 



These solutions are phase-locked in the sense that the dif- 
ference in phases at two fixed locations in space does not 
change with time. They describe both the synchronous 
(il 7^ 0, fc = 0) as well as the traveling wave solutions 
(ri 7^ 0, A: 7^ 0). The traveling wave solution are thus 
phase-locked solutions in which each oscillator has the 
same frequency Vt but the phase varies monotonically 
along the ring. Note that k is an integer due to peri- 
odic boundary conditions, and the value of (J)q can be 
taken to be zero by a translation. Substituting Eq. (Is]) 
into Eq. (IT]) gives the following dispersion relation, 



rj = cj— 



/ G{z) sin \VtT„ 



\z\ + Trfcz] dz. 



e'^K^Tjn ( — fc^TT^ 



mj 



(-1 + e«) (k2 + {kiT - f7r,„)2) (ft;2 _^ (fcjr + f7r„)2) 

(-1)*^^ (-f7r„, (-fc27r2 + k2 + f72^2j Cos[f7r„] - n (k^n^ + k"^ + n'^T^) Sin[f}T™]) 
(-1 + e«) (k2 + (ki: - riT„,)2) (k2 + ^kn + r!r„02) 



(6) 



r 



The solutions of the above general dispersion relation rep- 
resent the frequencies of various collective states of the 
system and its transcendental nature implies that fl can 
be multi- valued in principle for a given set of parameters 
a;,T,„, k and k. It is interesting to take various limits of 
Eq. (pi) and relate them to some of the past work on sim- 
plified models. In the limit when k — > 0, Eq. ^ reduces 
to, 



Q = LU — 



(1 - (-l)'=cos(r!T„))Or„ 



(7) 



-(fc7r)2 + (r!T™)2 
This is the reduced model that Zanette |^ had investi 



gated to obtain propagating waves in a system of globally 
coupled oscillators through the introduction of distance 
dependent delays. The other interesting scenario is that 
of K ^ 1, which corresponds to the situation of local cou- 
pling among the oscillators. In this case Eq. ^ reduces 
to. 



n 



-kn + Qt„ 



k-K + VLt„ 



+ (fcTT - Or„i)2 ^2 _4_ (J^^ ^ Q.Tm)'" 



(8) 

The local limit can also be approached by taking L ^ oo 
in the unsealed form of the model and hence of Eq. ^ 
in which r^ is replaced by L/v and k by kL. The mean 



delay (too) in this system equals -^ and Eq. Mm can be 



rewritten in terms of Too as 



O = w — 



1 



-p + V,To. 



+ 



P+ f^To, 



2 Vl + (P - ^^"00)2 1 + (P + f^Too) 



(9) 



where p is a real number denoting the wave number for 
the infinite system. This is a third-order polynomial 
equation in il and can at most have three real solutions 
in contrast to the higher number of multiple roots of the 
transcendental Eq. ^. Crook et al. [7 further simpli- 
fied the model by replacing the explicit delay time depen- 
dence by a space-dependent phase shift in the coupling 
term and thereby Eq. ^ further reduces to 



n = Lo — 



1 



2Vl + (p-roo)2 l + {p + f^y 



(10) 



In this approximation, fl is simply rescaled by the value 
of u! and was taken to be zero by Crook et al.[7] as w 
played no role in the stability of the system under this 
approximation. We now turn to the generalized disper- 
sion relation Eq. ([6| and discuss its solutions. 

The dispersion relation given by Eq. mf can be recast 
in the form : 



where 



H{nf,K) 



n 



L 



■ LO 



H{nf,K) 



(11) 



G{z) sin {ct,nf\z\+TTkz)dz (12) 



Fig. [1] plots the numerical solutions Q of Eq. (Ill as a 
function of f for k = 2.0 & w = 1.0 for different values 



of the wave number k. We note from Eq. (11) that the 



difference between the phase-locked frequency (51) and 
the intrinsic oscillator frequency (w) is a function of the 
corresponding Qf and the non-locality parameter k and 
has no explicit dependence on uj. This provides some 
simplification in the sense that if we plot the solutions 
of Eq. (11) in the phase space of fi — a; versus fif or 



equivalently in the phase space of H(flf, k) versus ilf, 
the results hold for any uj. Figure [2] plots these universal 
curves which are the solutions of Eq. ( |TT| ) for traveling 
waves with fc = 1, 2 and 3 respectively and for three dif- 
ferent values of k. For completeness we have also included 
the results for the synchronous state {k — 0) that were 
previously obtained in [15'. The three different values 
of K have been chosen to probe parametric domains in 
the vicinities of global coupling {k = 0.05, Fig. [2ja)), in- 
termediate (or non-local) coupling (k — 2.0, Fig. l2jb)) 
and local coupling (k — 10.0, Fig. [2jc)). While the equi- 
librium curves are universal (i.e. independent of uj) the 
stability domains of these solutions indicated by the solid 
portions of the curves do depend on uj. In Fig. [2]the solid 
portions of the curves in all the three panels denote sta- 
ble states for w = 1. The insets in each panel indicate 
the widths of the stable regions in ilf space for the cor- 
responding k values. In the next section, we discuss the 
method by which we determine these stable regions. 
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FIG. 1: (Color online) Synchronous frequency Q. for various 
wave mode numbers k for fixed values of k and oj. fc = 0, 1, 2 
and 3 are shown in blue, red, green and in magenta respec- 
tively and are also marked with the corresponding k values. 



III. STABILITY OF THE PHASE-LOCKED 

STATES 

The linear stability of the solutions 0si.fe(a;, t) of Eq. (fTl) 
is determined by the variational equation 



d_ 
dt 



u{x,t) 



/ G{z)cos[c,iflT\z\ + Trkz] 



where u{x,t) 

u[x,t) 

value equation 



X [u{x, t) — u{x — z,t — \z\cKT)]dz (13) 

= 4>{x,t) — (pn.ki^jt)- With the ansatz 
gAtgjTrnx^ A e C, 71 e Z, we obtain the eigen- 



/(A) = A + / G{z) cos {c^nf\z\ + nkz) 



j_ _ g-A|^|c„r^ 



dz = 



(14) 



Writing A = A^ -1- iXj and separating Eq. ( 14 1 into its 
real and imaginary parts we get 



X, 



-1 
1-e 

1 



X e 



G{z) cos (c„51r|z| -I- nkz) 

\H\z\c^f ^,Qg (Xj\z\c^f + nnz)] dz, (15) 

G{z) cos (cK51-f|2;| -I- nkz) 

^ sin{Xj\z\ci^f + nnz) dz. (16) 



-1 

-A_n|z|c„ 



The linear stability of the traveling wave state requires 
that all solutions of Eq. ( 14 ) have Aii; < for all non-zero 



integer values of n. The implicit and transcendental na- 
ture of the eigenvalue equation makes it difficult to obtain 
an analytic solution so we adopt a multi-pronged numer- 
ical approach to determine the marginal stability curves. 
For a given phase-locked state (i.e. fixing a value of k 
to 1, 2,3 etc) and a fixed intrinsic oscillator frequency 



K=0.05 




FIG. 2: (Color online) Plots of SI — oj vs Qf for various phase- 
locked states (fe — 0,1,2,3) at three different values of k. 
The shapes of all the curves are invariant to uj but the stable 
region is uu dependent. The solid black segments denote stable 
states for cj = 1. The thin segments shown in red denote the 
unstable regions. The insets highlight the width of stable 
regimes for different k values. 



w, we solve Eq. (11) for fl for a range of k values to 



span through the whole range of global to local coupling. 
We vary k from 10^^ to 15. The solutions so obtained 
are tested for eigenvalues with positive real parts by the 
Cauchy's argument principle which states that, the num- 
ber of unstable roots m of /(A) is given by : 



1 

27ri 



/(A) 



dX 



where the closed contour C encloses a domain in the right 
half of the complex A plane with the imaginary axis form- 
ing its left boundary. For each of the k values, r^ (and 



hence f) is varied from zero to the values where Xn = 
transitions are noted. Although this analysis is carried 
out for a large number of non-zero integer values of n, it 
is found that for synchronous solutions (i.e. k = 0), the 
lowest mode number, namely, n = 1 is the first one to 
get destabilized |15j whereas for the higher wave numbers 
{k = 1, 2, 3), the first mode number n to get destabilized 
is normally less than 10 for the parameter range covered 
in the present work. The analysis provides us r2,T,„ and 
n near A/j — transitions for any given values of n, uj 
and k. The transition values of fl and r^ are further 
refined (by using these values as initial guesses) by solv- 
ing simultaneously a set of three equations (Eqs |6|15 and 
16 1 with Afl — for fl,Tm and A/. We further find that 
A/ 7^ at the marginal point for any of the traveling 
wave solutions given by fc = 1, 2, 3.. indicating that they 
lose their stability through a Hopf bifurcation. This is 
in contrast to synchronous states (fc = 0) which were 
shown to lose their stability through a saddle-node bifur- 
cation with A/ = [TS]. One common feature that the 
traveling waves share with the synchronous states is that 
their stability regions are always restricted to the low- 
est branch of the equilibrium solutions where the curves 
have a negative slope. This is clearly seen in Fig. [2] and 
suggests that one can write down a heuristic necessary 
condition for the stability of the traveling wave solutions 
to be iJ' < 0, where the prime indicates a derivative of 
_ff (f2f, k) w.r.t Qf. Using the expression for H given in 
(12 1 the heuristic condition can be expressed as. 



\z\G{z) cos {c^nf\z\ + TTkz) dz > (18) 



For k = the above necessary condition applies to the 
synchronous states as was discussed in [15]. As was also 
shown in |15j , one can further obtain an analytic expres- 
sion for a sufficient condition for stability of the syn- 
chronous states given by, 



G{z) cos {c^nf\z\) [1 - cos (ttz)] dz > (19) 



This expression is easily obtained from Eq. 



(15 1 by ex- 
l)is 



ploiting the facts that the lowest mode number {n 
the first one to get destabilized and also the fact that the 
synchronous states lose their stability through a saddle- 
node bifurcation for which A/ = [T5] . 

The above expression provides simple and handy an- 
alytic stability criteria in limiting cases, e.g. fir < — ^ 

for global (k — > 0) coupling and fif < 4= for the lo- 
cal (k — > oo) coupling. In terms of frequency depres- 
sion (fi — cj), these inequalities for the stability become 



(17) n~Lu > - 



-"'(^) 



for global and H. 



UJ > 



73 



for 



2V2 4 

local coupling. We however note that under the assump- 
tion of constant delay t, the synchronization condition 
from Eq. ( 19 1 simply becomes fir < 7r/2 which agrees 



with the previously obtained results for constant delay 
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FIG. 3: (Color ) Stability diagrams in the parameter space 
of Sir vs hi for mode numbers k = 0,1,2 and 3, shown in 
blue, red, green and in magenta respectiveley. The solid black 
curves in each panel represent the condition H' — 0. Solid 
color areas represent stability regions for lj = 1 and the lighter 
shaded areas show the expansion of the stability region when 
cj — 10. Note that the stability region in panel (a) for A: = 
is independent of lv. 
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FIG. 4: (Color) (a) The superimposed stability diagrams of 
k — 1 and 2 for u = 1 shows a forbidden region where no 
phase-locked states are possible. The inset gives an enlarged 
view of the forbidden region. Similar forbidden regions are 
shown for uj = 0.8 in (b) & (c) and for lj = 0.5 in (d). The 
color coding is the same as in Fig. |3] 



systems [HI HI]- The stability condition obtained for lo- 
cal coupling (rif < 4=) matches with that of given by 
Crook et al. [7] for infinite system with non-local coupling 
incorporating propagational delays as phase shifts. 

For our present stability studies of traveling waves it 
is not possible to derive such an analytic expression for 
a sufficient condition since the transition to an instabil- 
ity is now through a Hopf bifurcation and hence A/ 7^ 0. 
Nevertheless the heuristic necessary condition ( jlS] ) serves 
a useful purpose by providing an upper and lower limit 
for the stability domains as we will see soon. Fig. [2] also 
shows some other interesting features, namely, the pres- 
ence of multi-stable regions where more than one trav- 
eling wave modes are stable (seen very prominently in 
(Fig. I2jc)) and the existence of a forbidden region (be- 
tween fe = 1 and fc = 2 for K = 2 in Fig. I2jb)) where no 
phase- locked states can exist. The size and shape of these 
regions are further found to depend on the magnitude of 
the intrinsic frequency w and of the non-locality parame- 
ter K. To investigate this dependence we have examined 
the stability of the traveling waves for two different val- 
ues of uj, namely, uj ^ 1 and a; = 10 and the results are 
displayed as a consolidated diagram in the phase space 
of rif versus k in Fig. [3] In this figure we have also plot- 
ted the H' = curves to demonstrate the fact that the 
stability regions in all cases are bounded by this heuristic 
necessary condition. In order to illustrate the sensitivity 
of the stability region to the value of the intrinsic fre- 
quency we have superimposed the results for w = 1 and 
a; = 10 in Fig. 3lb,c,d). The darker colored solid regions 
show the stable regions for uj — 1 and the superposed 
lighter shaded regions pertain to uj = 10. We note that 
a higher intrinsic frequency uj expands the stable region. 
The results in Fig. IS^a) do not depend upon uj. We also 
note that stable traveling waves are possible even with 
zero delay beyond some critical value of k which in turn 
depends on the particular wave mode number. This is 
clearly seen in Fig. l3|b,c) for fc = 1 and k — 2 respec- 
tively. 

The stability regimes of the k = 1 and k — 2 for cu = 1 
have been superposed in Fig. |4](a) highlighting the for- 
bidden region in the inset. In Fig (Hrb,c,d)) we have fur- 
ther explored the dependency on oj by superposing the 
stability diagrams obtained for uj = 0.8,0.8,0.5 respec- 
tively. An interesting development is the constriction of 
the k = 1 stability regime in the neck region seen in 
FigH[a) to its breakup into two disconnected regions (in 
Fig.TlJb)) to the total disappearance of the stable region 
in the low n regime in Fig. Illd). In Fig.llFd) we also no- 
tice disconnected stability regimes for k — 2 and k — 3. 
These regions exhibit a tongue structure for high con- 
nectivity (low k) whereas they submerge into the stable 
region of the synchronous state for low connectivity (high 
n). 

An interesting feature of the stability diagram is 
the existence of forbidden regions where both the syn- 
chronous as well the traveling wave states are unstable 
indicating a forbidden zone in parametric space where no 
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FIG. 5: Spatio-temporal cliaracteristics of a non-stationary 
breather state. Panel (a) shows the snap shot of the spa- 
tial pattern of phase after the transients are over. The 
corresponding inset shows the spatial variations of phase 
velocity(0) and the amplitude R of the complex order pa- 
rameter (Z). The panels (b) & (c) exhibit temporal behavior 
of R and (j) respectively. The parameter values are : u — 1, 
n — 2 and r™ = 6.4. 



stable phase-locked modes can be sustained. This is first 
shown in Fig. ^b) and subsequently in more details in 
Fig.ga-d). 
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FIG. 6: The same as in Fig. |5] but with parameter values: 
Lj = 1, K = 10 and Tm = 12.8. 



has a spatial structure and also displays regular periodic 
variations in time. Fig. pta) shows a snap shot of such a 
state with (j){x, t) as a function of x and the correspond- 
ing inset shows the spatial structures of (f>{x, t) and the 
amplitudes R of the complex order parameter Z defined 



Z{x,t) =Re"^ = I 



Q(^^ym^-^-t-\Ar^)\ rf^ (20) 



Z is a useful measure of the phase coherence of the oscil- 
lators in the system. Its amplitude < i? < 1 measures 
the system's coherence and $ is the average phase. 

Fig. p[b,c) show the periodic temporal behavior 
of R and (p respectively for this novel non-stationary 
breather state of the system. Fig. [6]^a,b,c) show another 
non-stationary breather state for different set of system 
parameters. Note that the spatial characteristics of 
the non-stationary state are intermediate between a 
traveling wave state (which would have a linear de- 
pendence of the phase on distance) to a chimera state 
(whose phase vs space curve has a nonlinear nature 
that is broken up with regions of incoherent regions). 
The non-stationary solutions have a smooth nonlinear 
curve for the spatial dependence of the phase (as shown 
in Fig. 5(a) and Fig. 6(a)). Its phase velocity (at any 
given spatial location) is a periodic function in time 
(as shown in Fig. 5(c) and 6(c)) in contrast to the 
traveling wave which has a constant phase velocity and 
the chimera state whose coherent portion also has a 
constant phase velocity. We believe these non-stationary 
breather states to be the intermediate link between 
the stationary traveling wave states and the clustered 
chimera states observed in |14) . in close analogy to the 
evolutionary sequence observed in no-delay systems 
between synchronous states, breathers and chimeras 
[12] . Indeed as we move away from the forbidden regions 
into the stable regimes for the traveling wave states we 
also find the co-existence of clustered chimera states in 
agreement with previous findings in [14j . 



CONCLUSIONS 



IV. NUMERICAL SIMULATION OF THE 
MODEL 



As an additional check of our stability results we have 
also carried out detailed and extensive numerical simula- 
tions of the discretized version of Eq. (IT]) and confirmed 
the existence of the stable and unstable regions shown in 
the various figures. In the course of these simulations we 
have also discovered the existence of a new class of non- 
stationary solutions of the system that have not been 
observed before. For these states the phase (p{x,t) and 
its temporal derivative (f>{x,t) have a non-linear depen- 
dence on space and the order parameter of the system 



To conclude, we have carried out a detailed stabil- 
ity analysis of traveling wave solutions for a generalized 
Kuramoto model of coupled identical phase oscillators 
that includes a spatially decaying coupling kernel and a 
space dependent time delay parameter. The comprehen- 
sive stability diagram presents a rich mosaic of multi- 
stable regions interspersed with unstable forbidden re- 
gions in the parametric space of the normalized quanti- 
ties representing time delay f , non-locality scale length 
K and the intrinsic frequency of an individual oscillator 
uj. The present study complements our previous work 
on synchronous states to provide a consolidated stabil- 
ity diagram of phase locked states for the generalized 



Kuramoto model. It also ties up some of earlier equilib- 
rium and stability results obtained from reduced models 
or simple limits of the present model to provide a more 
complete picture. Our investigations have also brought 
out important new findings such as the existence of for- 
bidden regions and the possibility of exciting stable trav- 
eling waves even with zero delay beyond some critical 



value of K. Likewise the non-stationary breather states 
are a novel find that provides further interesting informa- 
tion about the collective dynamics of this system and in 
conjunction with our stability results on traveling waves 
and synchronous states may prove useful in the predic- 
tion and interpretation of pattern formation phenomena 
in diverse practical applications. 
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